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Abstract 

The angular distributions of photoelectrons from atomic oxygen, 
nitrogen, and carbon are calculated. Both Hartree-Fock and Hartree-Slater 
v Herman~Ski liman) wave functions are used for oxygen and the agreement is, 

i 

excellent; thus only Hartree-Slater functions are used for carbon and 
nitrogen. The pitch angle distribution of photoelectrons is discussed and 
it is shown that previous approximations of energy independent isotropic 

. 2 

or sin 9 distributions are at odds with our results, which vary with 
energy. This variation with energy is discussed as is the reliability 
of these calculations. 
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I. Introduction 


the process of photoionization plays a significant role in many 
aeronomic and astrophysical phenomena [Dalgarno, 1967; Takayanagi and 
Itikawa, 1970]. If the process occurs in a dense medium, the only important 
quantity is the photoabsorption cross section; the angular distribution of 
the photoelectrons is inconsequential since it rapidly becomes isotropic 
due to collisions with the atoms and molecules of the medium. In a less 
dense medium, where the photoelectrons can travel considerable distances 
before losing their energy, the initial photoelectron angular distribution 
can be important. An example of such a situation occurs in the upper F2 
layer of the earth's upper atmosphere where photoelectrons produced by 
solar radiation ionizing atomic oxygen (produced by solar photodissociation 
of Og) can either deposit their energy locally or spiral along the lines 
of the earth's magnetic field and deposit their energy nonlocally at the 
conjugate point [Mariani, 1964; Whitten, 1968; Fontheim, et aU 1968; Kwei 
and Nesbit, 1968], i.e., magnetic focusing of the photoelectrons by the j 
earth's magnetic field. 

In this paper, the angular distribution of photoelectrons from atomic 
oxygen is investigated using Hartree-Fock (HF) wave functions. A previous 
work [Kennedy and Manson, 1972a] used Hartree-Slater (HS) wave functions in 
an incorrect theoretical formulation. Here the correct formulation is 
used to compare HS. and HF results. Agreement between these results is so 
good that we have extended the HS calculations to atomic nitrogen and 
carbon as well. 
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In Section II a brief discussion of the theory of photoionization and 

i 

photoelectron angular distributions and of the method of calculation is 
given. A derivation of the formulae relevant to the conjugate point pheno- 
menon is presented in Section III. In Section IV, we present and discuss 
the results of our calculations while Section V contains some final remarks 
and recapitulates the major results of this paper. 



II. Theory of Photoionization and Photoelectron Angular Distributions 

> : 

If we assume the wave functions of an atomic system can be well approx - 

J 

imated by anti symmetrized products of single particle wave functipns 

i 

# j 

(labeled and for initial and final states respectively), jthen the 

general expression for photoionization of an electron from an (n£[)^ ^S+l^ 

state to an [{(njO q-1 2Sc+1 L c >, (e£')3 2S+1 L‘ state is given by [bates, 

1946] 
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Here I is the.ionization potential of the • L state of the atom 
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to the 2Sc+1 L r state of the residual ion core, e is the photoelectron energy, 

! 

R the Rydberg energy, a the fine structure constant, a Q the Bohr radius, 
the statistical weight of the. initial state, t' the final photbelectron 
angular momentum, z > is the greater of £ and £', and c the relative multi pi et 
strength. The overlap integral y and the radial dipole matrix element 
R A * ( e) are given by 


Y = n 
passive 
electrons 




P n*W p n>;) d r | 2 


( 2 ) 


and 


«»■(*> 


p ^(r)r p ^,(r)dr 


(3) 


3 



where the continuum wave function is normalized such. that 

P^i(r) -* 7T - ^e” ^ sin(e^r + E~%ln2E%r - hli\. + o^' + 6^») s (4) 

r -+■ “ 

where o^« = argrU* + 1 - ie - ^) is the Coulomb phase shift and where 5^» 
is the phase shift with respect to Coulomb waves. This normalization of 
P^j,i is the usual normalization of continuum wave functions per unit 
energy range. 

The single particle wave functions for the initial discrete states 
were chosen to be Hartree-Fock (HF) functions compiled by Clementi (1965); 
the single particle discrete functions for the positive ion cores were 
similarly chosen to be the HF functions of Clementi (1965). The wave 
function for the final continuum electron is obtained in the field of the 
HF positive ion core. The numeri cal details are given by Manson and Cooper 
(1968) and Kennedy and Manson (1972b) and the continuum HF equations have 
been presented in detail by Dalgarno, Henry and Stewart (1964). 

In addition Herman-Ski liman . (HS) wave functions have also been employed 
in these calculations for comparison with the HF results. In using these 
wave functions no core relaxation is considered and the continuum function 
is found from the central potential appropriate to the ground state of the 
atom. A full discussion of this type of calculation is given by Manson and 
Cooper (1968). 

9C+] 

The angular distribution of photoelectrons from an initial L state 

pc +1 

to a c L c state of the residual ion core is given for unpolarized light 
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by [Cooper and Zare, 1968, 1969] 

/ 

dtf-j(e) cr-j ( e) B^e) ( 5 ) 

= -47" [1 ■ T - p 2 ( cose )] 
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where j 

1 

A+l | 

° i(e)= L‘ £ '= £ - 1 a n* (LS ’ L c S c> cA ,L< ) ( 6 ) i 

i 

0 is the angle between the incident photon beam and' photoelectron direction' 

2 

P 2 ( x ) = %(3x -1), and 3-j(e) is the asymmetry parameter. The form of the 

| 

angular distribution [Eq. (5)] can be obtained from quite general consider- j 
ations depending only on the absorption of radiation by an unpolarized target oc- 
curring via an electric dipole process [Yang, 1948]; deviation from the form of this 
/ equation implies the presence of absorption via processes other than electric dipole. 
The asymmetry parameter, ^-(e), however, does depend upon the details of the 
calculation. Specifically, for the LS coupled anti -symmetric products of ^ 
single particle Hartree-Fock functions we are using, the asymmetry parameter 
is a sum over the s's corresponding %o each of the possible values of the 
angular momentum j t transferred to the atom by the photon, each of these 
B's being weighted by the relative cross section for each value of j t 
[Dill and Fano, 1972; Fano and Dill, 1972; Dili, 1973] 

6 i = ? 3(j t )o(j t )/z c( j t ) . ( 8 ) 

J t h 
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Transitions where the parity change of the target is (-l) Jt are said to be 
parity favored and 3(j t ) is given by [Dill and Fano, 1972; Dill, 1973] 




t'fav 


(j t +2) |S + (j t )| 2 +(j t -l)|Sjj t )| 2 -3[j t (j t +l)]^[S + (j t )^(j t )+S*(j t )S.(j t )] 


(2j t+ l)[|S + (j t )| 2 + |S_(j t )| 2 ] 


(9) 


where S ± (j t ) denote the photoionization amplitude for a given and for 
i' = j t ±l; values of for which the parity change is -{-l)^t are called 
parity unfavored and have a sin e distribution [Dill and Fano, 1972; Dill, 
1973], i .e. , 


^t^unf “ » 
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independently of dynamics. The integrated cross sections a(j t ) have the 
structure [Dill, 1973] < 
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where S 0 (j t ) is the photoionization amplitude for the value of 


6 


and * is the photon wavelength (divided by 2n). Detailed discussion of 

/ \ 

these matrix elements is given by Dill, Hanson, and Starace (1974). As 
an example of the results, 3 -j for the P -*• D photoionization (labeled 
by i) for oxygen is given by 


3 = 


_ 3|?j |i (l)| 2 -3^[5'i i 1 (l)I- ii (l)*'fcc]-5[|S^ i (2)|2-|? d|i (3)|2] 
3[|^ >i (D| Z +|Si ii (l)| 2 ]+5|S^ 1 -(2)|2+7|Si ji (3 )| 2 

(13) 

where cc is a complex congugate and 
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(14c)* 


P. 3, 
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in which a s and are the coulomb phase shifts defined previously. The 
superscripts on the phase shifts and dipole matrix elements refer to the 
different total angular momentum states of the photoelectron-ion system 
after photoionization: e.g., a d-wave photoelectron plus a ion core can 
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couple to possible final states 3 $, 3 P, and 3 D. 

Note that the second term in the numerator of Eq. (13) introduces a 
cosine of the difference of the phase shifts between the d- and s-continua 
very much like the simple Cooper-Zare expression [Cooper and Zare, 1968, 
1969]. The other terms, which do not appear in the Cooper-Zare formulation, 
introduce the cosine of phase shift differences between the various possible 
d-continua (^S, ^P, or ^D). However, it is seen from Eqs. (12b - d) that 
if the phase shifts and dipole matrix elements of d-continua are the same, 
(as they are in a HS type calculation), the j t = 2 and j t = 3 contributions 
vanish and Eq. (11) reduces to the Cooper-Zare formula for an initial 
p-state. 


0 


i 



R d - 2R s R d cos Ud - e s > 



+ 
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where ^ + i = + ] + ± i • Eq (15) shows that within the HS approxi-j 

mation, 3-j (e ) is the same for each channel and depends only upon the 
photoelectron energy, s. Thus when a HF calculation results in differing 
matrix elements and phase shifts for differing final state multiplets with 
the same ion core, the Cooper-Zare formula will be in error to an extent 
determined by the magnitude of these phase shift differences. In closing 
this section we emphasize that Eq. (15) applied to the angular distribution 
of ill photo-electrons in a transition from a given atomic state to a given 
ionic state. The discussion of this point in a previous paper [Kennedy and 
Manson, 1972] is incorrect and should be disregarded. The a results in 
that paper for the ->■ ^P channel is the correct one for all channels in 
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the Cooper-Zare approximation, whose validity for oxygen is examined below. 
> '■ 

III'. Pitch Angle Distribution of Photoelectrons 

For the conjugate point phenomenon, we are interested in the angular 

distribution of photoelectrons with respect to the earth's magnetic field 

direction, i.e., the pitch angle distribution. To get this, we have the 

situation shown in Fig. 1 where the magnetic field is in the z-direction 

f 

and the incident photon is in the x-z plane, making an angle y with the 

/ 

z-axis. Assume the photoionization takes place on the z-axis with the 
photoelectron going in the direction shown. The angular distribution of 
photoelectrons is given generally by 


da 4 tt 

f(cose) = s A n P n (cose) 7 A n z Y™(y,0) Yf}(e',$) 

n n 2n + 1 m = -n 




(16) 


using the addition theorem for spherical haronmics. We are interested in 

‘.the pitch angle distribution, ^7 7 without regard to <f>. Thus, integrating 

d(cose') 

over $ we obtain 


da 


d(cose') n 


= 2tt z A n P n (cos’r) P„(cose* ) . 
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For an electric dipole interaction Ag = o/(4tt) and A 2 = -(b/2)Aq so that 


da 

d(cose‘ ) 


a 

- [1 - I P 2 (cos Y ) P 2 (cose' )] 


which is the pitch angle distribution. 


( 18 ) 
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IV. Results and Discussion 


> i 

Calculations have been performed for photoionization of the 2p subshell 

3 1 1 

of oxygen from the ground configuration in the P, S, and D states in both 
HF and HS approximations as discussed in Section II. The calculated cross 
section for photoionization from the 3 P state of the ground (2p)^ configuration 
of neutral oxygen is given in Fig. 2. The results show fairly good agreement 
among the HF-length (HFL), HF-velocity (HFV), and HS formulations. In addition, 
our results agree essentially exactly with the HF results of Dalgarno et al 
(1964). This is, of course, to be expected and serves to confirm our compu- 
tational technique. 

More importantly, it has been shown that the angular distribution of 
photoelectrons is only very weakly a function of the cross section; it is 
the phase shifts of the continuum waves which are all important [Manson and 
Kennedy, 1970; Kennedy and Manson, 1972b]. This is true in any spectral 
region where neither a shape resonance or a Cooper minimum occurs; we note 
that neither is in evidence for the photoionization of oxygen. Further, 
the phase shifts of the various ed's in the HF and HS approximation are 
essentially the same. Thus it is expected that the angular distribution 
parameters, the 3-j f s , will not differ very much in the various approximations. 
This is borne out by Fig. 3, which shows the b's for the various possible 
photoionization transitions computed in the HFL and HS approximations. The 
agreement among the several approximations is seen to be excellent. Thus 
we conclude that for light atoms, the HS approximation wil predict excellent 
B's and the Cooper- Zare formulation is adequate. 

In previous work relating to the conjugate point phenomenon, the angular 
distribution parameter was assumed to be isotropic (3 = 0) by Fontheim et al 
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(1968) or sin^e (3 = 2) by Mariam' (1964) and Whitten (1968). From Fig. 3 
however, it is apparent that the 6‘s vary with energy. Further, an isotropic 
distribution is valid only just above threshold, and 3-j never gets near 2. 

For energies a bit above threshold and thereon, 3-j ^ .1 or about halfway 
between isotropic and sin^e. Hence it appears that the approximations made 
in connection with the angular distribution of photoelectrons from atomic 
oxygen are not valid in any spectral region, and even the assumption of 
constant 6 is unwarranted, based upon our results. 

The variation of 3 with photoelectron energy s is 
caused principally by the variation of the phase shift difference 5 d - £ s > 
with energy. This difference is made up of the sum of the Coulomb phase 
shift difference given by [Manson, 1973] 

1 1 

a d " °s = 2/e -tan“ /e (19) 

and the non-Coulomb 6 d - S s . For oxygen at threshold <5 d * tt, while from F.q. 

(19) a d - a s = -tt, so that the phase shift difference at threshold is j 

about -2tt. Thus the cos term in Eq. (15) is about unity so the two terms •’ 

in the numerator tend to cancel, leading to the near zero threshold value 

{ 

of 6 shown in Fig. 3. As the energy increases to e = 1 Ryd. the only 

/ 

major change is in the Coulomb phase shift difference which becomes 'v -7r/2 
so that the C d ' “ _37T / 2 and cos term in Ec l- 0 5 ) vanishes; 3 is 

therefore about unity. This rapid rise of 3 just above threshold is a 
general characteristic of atomic photoelectron angular distribution 
parameters and is caused by the Coulomb phase shift variations in this 
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energy region [Kennedy and Manson, '1972b; Manson, 1973]. At higher energies, 
no major variation in matrix elements or phase shifts occur, so 3 varies 
only very smoothly as shown. 

. The adequacy of the HS approximation, as discussed above, for low-Z 
elements gives impetus to performing further calculations using this 
approximation. To this end, then, we have calculated the angular distribu- 
tions of the photoelectrons from the states arising from the ground config- 
urations of carbon and nitrogen. 

The 3 1 s for N and C are shown in Fig. 4. For N, the results are 
substantially like 0. The 3-j for N is substantially like that for 0 since 
the dipole matrix elements and continuum wave phase shifts are fairly close. 
The variations of 3 -j with e, . therefore , are results of effects discussed 
in connection with oxygen previously. 

For carbon the results show effects substantially the same as in 0 
and N and for exactly the same reasons as discussed above. It is thus clear 
that the approximation of a constant g = 0 or a = 2 is quite poor in these 


cases as wel 1 . 


V. Final Remarks I 

i 

From the previous section it is seen that the photoelectrons from atomic 
oxygen have an initial angular distribution which is neither isotropic nor 
sin 0 . Further, we find the initial angular distribution to be strongly 
dependent upon the energy (wavelength) of the impinging radiation. Since the 
distance travelled to the conjugate point is larger than the mean free path 

of the photoelectrons, it was assumed by Fontheim et aj_ (1958) that the photo- 

/ 

electron angular distribution is rapidly isotropized. This assumption, however, 
is invalidated by the fact that the electron collisions are primarily elastic 
which are generally sharply peaked in the forward direction [Mott and Massey, 

1965]. Thus, although collisions do occur, most of them do not affect the 
angular distribution significantly. We conclude, then, that the initial 
photoelectron angular distribution is significant for the conjugate point 
phenomena. 

As to the accuracy of our calculations, comparison with experiment shows 

excellent agreement [Kennedy and Manson, 1972; Mitchell and Codling, 1972; 

Lynch, Gardner, and Codling, 1972; van der Wiel and Brion, 1973] for noble'* 

gas atoms. Some modification of the theory is necessary when fine structure 

\ 

is significant in the initial state of the atom as pointed out by Dill (1973), 
but these effects are expected to be inconsequential here. This is presently 
under investigation. In addition, this work does not consider the effects of 
autoionization which will affect the angular distribution in the vicinity of 
an inner shell threshold [Dill, 1973], in particular the 2s. This affects things in 
only a very small energy range, a few eV below the 2s threshold, sc it should not 
be too important for the conjugate point phenomenon. Finally, it has been 
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shown theoretically by Amusia et aj_ (1972) that the inclusion of correlation 
in initial and final states significantly affects angular distributions only 
near Cooper minima [Cooper, 1962] and for subshells whose cross section is 
strongly perturbed by a nearby subshell whose cross section is very much 
larger via the interchannel interaction [Fano, 1961]. It is thus felt that 
the theoretical B's presented herein are of sufficient accuracy for applica- 
tion in aeronomic problems. 
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Schematic representation of the photoionization process. 

Photoionization cross section of 0( 3 P) -> o + ( 4 S 2 n 2 pl + p 
l he ^ S A lld , c ^ rve is the Hartree-Slater (HS) result ’and the ‘ 
vel length ’(HFL) and 

24fS{ c „™, 

channels'from SS 9 ^-'4 H ^.^? u !^.f or al 1 the Photo^onization 
configuration. 
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for 2p Photoionization of atomic 
N and C calculated using HS y/ave functions. 
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